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ABSTRACT 

We present very high resolution hydrodynamical simulations of accretion discs in black hole 
X-ray binaries accreting near the Eddington limit. The results show that mass loss, irradiation 
and tidal interactions all have a profound effect on the observed behaviour of long period 
X-ray transients. In particular, the interplay of all of these effects in the outer regions of the 
accretion disc is able to drive long time-scale (weeks to years) variability is these objects, and 
is a possible origin for some of the extreme variability of GRS 1915+105. 

Key words: accretion, accretion discs - binaries: close - stars: individual: GRS 1915+105. 



1 INTRODUCTION 

The outbursts of long period X-ray binaries involve accretion rates 
close to the Eddington limit. Consequently, the evolution of the 
massive accretion discs in these objects is expected to be highly 
influenced by mass loss and irradiation by the central X-ray source. 
Outbursts of this nature are thought to power the most luminous 
stellar X-ray sources such as the galactic microquasars and at least 
some of the recently discovered ultra-luminous X-ray sources in 
nearby galaxies l Kins 2002). In this paper we investigate the evo- 
lution of accretion discs and the origin of the long-term variability 
in the outbursts of long period X-ray transients, using the galactic 
microquasar GRS 1915+105 as an example system. 

The galactic microquasar GRS 1915+105 (V1487 
Aql) was discovered as an X -ray transient in 1992 
ICastro-Tira do. Brand t & Lundtn l992). and has been observed to 
be extremely luminous ever since. This binary system contains a 14 
Mq black hole I Greiner, Cuby & McCaushrean 2001) accreting 
from a late-type giant of mass 0.8 ± 0.5Mq ( Harlaf tis & Greinerl 
2004) via Roche lobe overflow. GRS 1915+105 is unique among 
accreting Galactic black holes spending mu ch of its time a t super - 
Eddington luminosities iDone. Wardzinski & Gierlin ski 2004). 
It is an extremely variable source, exhibiting dramatic, aperiodic 
variability on a wide range of timescales, from milliseconds to 
months. FigureQshows the Rossi X-ray Timing Explorer All-Sky 
Monitor observations of X-ray flux from GRS 1915+105 over the 
past decade. 

The fundamental origin of this spectacular behaviour is 
thought to be an outburst similar to those exhibited by other 
X-ray transients, but of a much longer duration. The outbursts 
of X-ray transients are believed to be initiated by a thermal- 
viscous accretion disc instability, and prolonged by the heating ef- 
fect of the central, irra diating X -ray sourc e l Kins & Ritten ll998t 
iDubus. Hameurv &La sota 2001 ; Truss et al. 2002). The extremely 
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long time-scale of the outburst of GRS 1915+105 is a consequence 
of its long 33.5 day orbital period, which implies a very large disc 
size (7?disc ~ 3x 10 12 cm). Such a large disc can accumulate 
an enormous amount of mass during quiescence, which is able to 
support a subsequent super-Eddington outburst over many decades. 

So far, most models of the variability of GRS 1915+105 have 
concentrated on explaining the origin of the short time-scale lumi- 
nosity changes associated with the inner regions of the accretion 
disc (for example, Belloni et al. 1997), although Kins et al. 1 2004) 
have proposed a mechanism for generating variability on a wider 
range of time-scales, in which outflows of mass from the disc are 
driven by localised magnetic dynamo processes . However, the ori- 
gin of longer-term, aperiodic variations on timescales of hundreds 
of days remain poorly understood as these are most likely driven 
by dynamical evolution of the outer disc. 

In this paper we present the results of two-dimensional (2D) 
calculations of the long-term evolution of an accretion disc in a 
long-period X-ray binary accreting close to the Eddington limit, 
using GRS 1915+105 as an example system. We include the ef- 
fects of tidal instability, central irradiation and mass loss on the disc 
evolution. We show that the interplay of these effects in the outer 
regions of the disc can induce similar variability to that observed in 
GRS 1915+105. 



2 NUMERICAL MODEL 

In this section we discuss several key refinements that have been 
made to an existing smooth particle hydrodynamics (SPH) code. 

SPH is a Lagrangian scheme used to describe the dynamics of 
fluid flow. A continuous gas is described by an ensemble of par- 
ticles moving with the local fluid velocity. Local properties such 
as velocity, density and temperature are determined by an inter- 
polation over the neighbouring particles and in this way the fluid 
equations may be solved. In the continuum limit (that is, for a large 
enough number of neighbour particles), the SPH equations for mo- 
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Figure 1. «XT£ All-Sky Monitor (ASM) X-ray flux observations of 
GRS 1 9 1 5 + 1 05 . The system shows a wealth of variability over a wide range 
of time-scales and has been in a continuous state of outburst since 1992. 
RXTE results are provided by the ASMJRXTE teams at MIT and at the RXTE 
SOF and GOF at the NASA Goddard Space Flight Centre. 



mentum, thermal energy and so on reduce to the exact fluid equa- 
tions. The reader is referred to MonaahanJ I1992T) for a review. 

The code used for this work is based on a vectorised version 
of the three-dimensional accretion disc code developed by Murrav 
(1996). There are four key elements of physics that must be in- 
cluded in a code that attempts to describe the behaviour of a sys- 
tem like GRS 1915+105, accreting near the Eddington limit: an 
implementation of the viscous instability to drive the outburst, a 
description of the tidal forces of the secondary star on the disc, the 
effect of self-irradiation of the disc by the X-rays generated near 
the primary and mass loss when the Eddington limit is breached. 
We will describe each of these in turn. 



sound speed (but allow a continuum of states between). However, 
as we shall see, once the outburst is underway the disc as a whole 
spends little or none of its time in a single state and most of its 
behaviour is governed by local processes; changes in sound speed 
and mass loss occur in relatively small regions, rather than globally 
across the disc. It is not clear, then, that a fully self-consistent calcu- 
lation - computationally, this remains inviable in a two-dimensional 
hydrodynamic simulation of this magnitude - would produce vastly 
different results. Until this becomes possible, we believe that the 
method we describe here, which allows for realistic local changes 
in sound speed, offers a better prospect than the fully-isothermal 
schemes of the past. 

The sound speed and temperature are changed in response 
to the local surface density calculated at the position of each 
particle. We define critical surface densities, E max and £ m i n 
to trigger an outburst and to return to quiescence respectively. 
These critical values scale almost linearly with radius in the 
disc I Ludwia. Mever-Hofm eister & RitteJ 19941) . When the surface 
density at a particle crosses one of these values, both the sound 
speed and the viscosity parameter are changed smoothly according 
to 

/,\ (Cs hot + Cs cold) i (Cs hot — Cs cold) , , / t \ ,~, 

Cs(t) = ± — '■ — '■ tanh I ttJ (2) 

and 
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The thermal time-scale for the change, £ t h, is defined as the ratio 
of the heat content per unit (projected) area of the disc, Ec^ , to the 
total rate of energy loss per unit area, which taking into account 
both faces of the disc is 9^£GMi/4_R 3 . Therefore, the thermal 
time-scale 
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2.1 Viscous instability 

The disc instability, which is driven by the onset of ionization in the 
temperature range 6,000 to 7,000 K, is implemented by a viscous 
switch based on the local disc conditio ns. Th e viscosity, v, is given 
by the formula of lShakura & Sunvaevl ll973l) : 



ac s H, 



(1) 



where c s is the local sound speed, H the scale height and a is the 
viscosity parameter. It is well known that the increase in temper- 
ature at the transition from quiescence to outburst is insufficient 
to provide the angular momentum transport required to drive the 
mass accretion rates observed during an outburst. Hence, it is in- 
ferred that the heating is accompanied by an increase in the viscos- 
ity parameter a. In previous SPH studies of outbursts in accre tion 
discs, a was allowed to change instantaneously l Murrav fll998t) . or 
on some fixed time-scale l Truss et al. 2000), but c s was kept fixed. 
Here, we allow the sound speed as well as the viscosity parameter 
to vary on a time-scale that is determined by the local conditions. 
Recently, a similar method was appl ied to the ou tbursts of cata- 
clysmic variables JTruss. Wvnn & Wh eatlev 2004). The method is 
able to reproduce a full, realistic outburst amplitude with none of 
the scaling factors required by the previous schemes that used a 
constant sound speed everywhere. In this approach, we would not 
reproduce the expected temperature profile of a disc in steady state, 
where c s oc R~ 3 ' s , because we have lower and upper bounds to the 



2.2 Tidal forces and mass transfer 

The calculations are performed in the full Roche potential of a bi- 
nary system. In this way the tidal influence of the secondary star is 
intrinsic to the scheme. For the mass ratio q — M2/M1 ~ 0.07 
inferred for GRS 1915+105, we expect the tides to have a signif- 
icant effect on the orbits near the outer edge of the accretion disc, 
and to launch a two-armed spiral wave. Mass transfer from the sec- 
ondary star is implemented by the injection of particles at the in- 
ner Lagrangian point. The rate of mass injection is independent of 
the physical conditions at the edge of the disc and is kept constant 
throughout the calculation, regardless of whether the disc is quies- 
cent or in outburst. 



2.3 Irradiation 

We use a simple model due to lKing & RitteJ Jl99ut) . in which the 
X-rays emitted near the black hole keep the disc in the hot, high 
viscosity state out to a certain radius Rh which is determined by 
the central accretion rate. Matching the irradiation temperature to 
the minimum temperature to keep hydrogen ionised, Th, we have 



-Rh 



rjMc 



47TO-T, 



T f(P,H,R) 



(5) 



where r\ is the accretion efficiency and / is a factor dependent on 
the albedo, the disc geometry and the form of the irradiating source. 
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For the typi cal parameters of a disc around a black hole given in 
Kins & Ritter 1 1998), we have 



R h = 1.9 x 10 



"(<&)' 
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(6) 



where Ms is the central accretion rate in units of 1CP 8 Mqvt -1 . 
In the calculations, all particles within Rh are given the high state 
viscosity Qhot and sound speed c s ,hot- 

2.4 Mass loss 

In the next section we will discuss the results of three calculations 
of outbursts in GRS1915+105. In each calculation, the mass loss 
due to the Eddington limit being exceeded is dealt with in a dif- 
ferent way. In the first, which is intended to act as a control case, 
the limit is ignored and Rh is allowed to evolve for accretion rates 
beyond the Eddington limit. In the second, Rh is only allowed to 
increase up to a maximum value corresponding to a central accre- 
tion rate equal to the Eddington rate. 

In the third calculation, individual particles are removed lo- 
cally from the accretion disc where the Eddington limit is breached, 
not just at the inner boundary. The local mass accretion rate at each 
radius in the disc is calculated from 



M(r) - -2irrv T (r)Y,(r). 
The Eddington accretion rate at a radius r is 
rL E dd 



M Edd (r) 



GMi 



(7) 



(8) 



where L^dd is the Eddington luminosity. Particles are removed 
from the disc before they cross the inner disc boundary wherever 
M(r) > M E dd(r),or 



M{r) > M E 



^ (i 



(9) 



where Ro is the radius of the inner boundary. 

It is of primary concern in such a scheme, where particles or 
groups of particles are being removed from the flow, that spatial 
resolution is not degraded severely. With this in mind, we use a 
smoothing length that varies according to the local density. The re- 
sponse of the smoothing length is such that each particle maintains 
at least 50 nearest neighbour particles. 



3 RESULTS 

In this section we present the results of three very large two- 
dimensional SPH calculations of accretion discs in a binary sys- 
tem with parameters similar to those of GRS 1915+105: P rb = 
33.5 days, Mi = 14 M Q and M 2 = IMq. We use Cs, co ld = 
0.05aJ! or b, C s ,hot = O.15af2 or b, which correspond to mid-plane 
temperatures of ~ 5000 K and ~ 45000 K respectively, assuming 
a mean molecular weight fj, = 0.6, appropriate for a fully ionised 
mixture of gases (strictly, of course, this value of fj, is not self- 
consistent for the neutral gas in the low state. A better estimate of 
fl in the low state would be higher, although this would have little 
effect on the global evolution of the disc through an outburst in the 
simulations). 

The accretion efficiency is set to 77 = 0.1. Our choice of vis- 
cosity parameters, a co id = 0.1 and Qhot = 1, is slightly higher 
than we would ideally choose, but is necessitated by our need to 
simulate a very long time series of data with a limited amount of 
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Figure 2. Top to bottom: energy dissipation rate, central accretion rate and 

irradiated radius as a fraction of binary separation during simulation 1 . 



dedicated supercomputer time. We discuss the implications of us- 
ing different viscosity parameters fully in Section 4. 

The mass transfer rate is calculated from the formula of lRitted 
(1999): 



M, 



M2 

M 



7.3 X10- 10 £i P o rb 98 (d) M yr-\ (10) 



which is appropriate for nuclear time-scale-driven mass transfer 
from an evolved companion star. In this case, the formula gives 
-M 2 ~ 2 x 10~ 8 MQyr~ 1 . This is already close to the Eddington 
rate (near the black hole); the implications of such a rate are also 
discussed in Section 4. 

Each of the three simulations have identical initial conditions: 
a steady-state accretion disc with iVo = 439, 603 particles. This 
places these simulations among the largest accretion disc simula- 
tions ever attempted with SPH. The simulations were run on the 
UK Astrophysical Fluids Facility (UKAFF) during allocated time 
awarded in late 2003. 



3.1 Simulation 1 

The accretion disc was allowed to evolve through an outburst with 
no restriction on the irradiation radius and no mass loss. Figure [2] 
shows the resulting time-variations in total rate of energy dissipa- 
tion per unit mass in the disc (summed over all particles), du/dt, 
central accretion rate and irradiated radius, 7?i rr . 

The central accretion rate rises by two orders of magnitude to 
apeak at 3.5 x 10 _7 MQyr _1 , and the entire disc is irradiated. The 
outburst profile that results is very simple: as the disc is kept in the 
hot, highly viscous state by the central X-rays, we see an exponen- 
tial decay as the classic 'flat-top' to the outburst. After six months 
or so, the irradiation radius has retreated to the inner edge of the 
disc, and we see a rapid decline to the quiescent state, which per- 
sists until the end of simulation. The number of particles remaining 
at the end of the outburst is 189,307, corresponding to the net accre- 
tion of 57% of the total mass of the disc. Each particle has a mass 
2.09 x 10 _13 Mq, so the total mass accreted during the outburst is 

5.2 x 1O" 8 M . 

Clearly, this simulation cannot explain GRS 1915+105, but 
the behaviour is consistent with the results of simulations of the 
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Figure 3. Top to bottom: energy dissipation rate, central accretion rate and 
irradiated radius as a fraction of binary separation during simulation 2. The 
brief gaps in the data in the middle panel correspond to times when the 
central accretion rate was zero. 



short period X-ray transient A0620-003 in iTruss et alj <2002l) . 
where a similarly simple outburst was obtained in the case where 
the entire disc becomes irradiated by the central X-ray flux. 

3.2 Simulation 2 

In this simulation, we assume that the Eddington limit holds at the 
inner boundary of the accretion disc. Material accreting at a rate 
higher than MEdd(-Ro) is assumed to be lost from the system and 
does not contribute to increasing the fraction of the disc that is illu- 
minated by the central X-rays. 

Therefore, the irradiated fraction of the disc is limited to a 
maximum value corresponding to accretion at the Eddington rate 
through the inner boundary. We use AfEdd(-Ro) = 10 _7 Mq yr _1 , 
which gives i?h,max = 0.15a, where a is the binary separation. 

The contrast between the resultant outburst, shown in Figure 
[5] with that of the previous simulation is stark. As expected with 
identical initial conditions, the initial rise to outburst maximum pro- 
ceeds exactly as before. However, after a short time, Rh reaches its 
maximum value, and remains there for about two orbits of the bi- 
nary while the central accretion rate is super-Eddington. Following 
this plateau, the accretion rate declines and enters a phase of aperi- 
odic variability that is reflected in the energy dissipation rate in the 
disc. We note from the lower two panels of Figure|3|that the maxi- 
mum central accretion rate during these flares does not exceed the 
Eddington rate at the inner boundary. The flaring behaviour con- 
tinues until the end of the simulation at t — 80 orbits, long af- 
ter the last data point shown in the figure. Of course, we do not 
know for sure if the system would return to quiescence if the sim- 
ulation could be run indefinitely. However, from the data we have, 
it seems that the system has found a quasi-steady state: the time- 
averaged accretion rate over the whole simulation is very close to 
2 x lO^MQyr" 1 , which is the rate at which mass is added to the 
disc from the secondary star. 

Figure|4| shows the evolution of the surface density profile of 
the disc during the initial stages of the simulation. While the cen- 
tral accretion rate is super-Eddington, a large spike of gas builds up 
at the hot/cold boundary. As the highly viscous, hot gas spreads, 




Figure 4. Azimuthally averaged surface density profile during the initial 
stages of simulation 2. The dotted lines are the critical surface densities 
£ m in and S max described in the text, although in all but the top and bot- 
tom panels everything inside R = 0.15a is kept in the hot state by the 
irradiation. The times are quoted in orbital periods. 



most of the gas outside R^ is pushed over the critical limit E max 
and starts to participate in the outburst. However, this plot does not 
seem to explain the origin of the flaring behaviour. The surface den- 
sity profiles are produced by averaging over azimuth, and much in- 
formation about the structure of the flow is lost in producing them. 
Figure [5] shows the full density structure of the flow at two times 
during the simulation; at maximum R^ at £=0.4, and later on during 
the flaring behaviour at £=5.0. The initial state of the disc is rather 
smooth with some weak spiral arm features (these can still be seen 
in the outer part of the disc in the left-hand image of Figure|5J. Af- 
ter the onset of the outburst, however, it can be seen immediately 
that the flow through the fixed hot/cold boundary is neither steady 
nor smooth. The disc is left with a clumpy, inhomogeneous struc- 
ture, and it is the accretion of the over-dense clumps of gas that lead 
to the subsequent flaring behaviour. This is a direct consequence of 
fixing i?h for a short time when the accretion rate onto the black 
hole is super-Eddington. We have already seen in simulation 1 that 
this structure does not develop if Rh is not fixed. 

It is interesting to note that the tidally-induced features are 
preserved as the outer edge of the disc expands during the outburst. 
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Figure 5. Surface density maps of the accretion disc in simulation 2 at £=0.4 orbits (left) and £=5 orbits (right). The colour scale runs from black at gem -2 
through red to yellow at 12 gcm~ 2 . The secondary star to the right of the mass stream and the inner part of the accretion disc(white) are not modelled. The 
tidally-induced spiral structure and non-circular outer orbits are preserved throughout the simulation despite the inhomogeneity of the flow. 



Despite the unevenness of the flow, a strong two-armed spiral wave 
is evident by £=5.0 orbits (right hand side of Figure|5} that pene- 
trates deep into the disc. In fact, it is the interaction of the spiral 
arms with the hot/cold boundary that leads to the variability in this 
simulation. While the gross structure of the spiral arms is preserved, 
locally their structure is determined by the fluctuations in pressure 
and density induced at Rh. 

The gas near the outer edge moves on eccentric orbits, and this 
eccentric structure does indeed precess in a prograde fashion in ex- 
actly the manner expected for a tidally-unstable disc. There is no 
evidence of superhumps: the periodic variations in dissipation rate 
observed in tidally-unstable cataclysmic variables that are caused 
by gas in orbits resonant with the orbital frequency of the binary. 
However, this is unsurprising because here the dominant compo- 
nent of the energy dissipation rate comes from the irradiated inner 
disc, not the tidally stressed outer orbits. 

The local mass accretion rate through the disc during the out- 
burst is extremely high. Figure |6| shows the evolution of the local 
accretion rate profile for the same time slices as the surface density 
plots in Figure|4| The initial steady-state profile is given in the top 
panel and is reproduced as a dotted line in the lower panels to aid 
comparison. The mass transfer rate at the outer edge of the disc re- 
mains constant at all times and save for some structure associated 
with the spiral arms near R=0.4a, the steady-state profile is con- 
stant with radius at a value close to 2 x 10 -8 Mq yr _1 - the rate 
of mass transferred from the donor star. When the spike of gas de- 
velops at the irradiated boundary, the local accretion rate rises to a 
maximum near 10 _6 Mq yr~\ well sufficient to sustain a super- 
Eddington central rate. By £=5 orbits, when the main part of the 
outburst is over and the flaring behaviour is underway, the profile 
approaches that of the steady initial state, but varies in response to 
the local inhomogeneity of the flow caused by the clumps and spiral 
structure. 



3.3 Simulation 3 

Here, mass is removed from the disc locally whenever the accretion 
rate exceeds the local Eddington accretion rate. This is the condi- 
tion that has been defined in equation|9| where we use Rq = 0.05a 
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Figure 6. Azimuthally averaged local mass accretion rate during the initial 
stages of simulation 2. The initial profile is reproduced as a dotted line in 
each case to aid comparison. The times are identical to those in Figured 
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Figure 7. Top to bottom: energy dissipation rate, central accretion rate and 

mass-loss rate during simulation 3. 



and A/Edd(^o) = 10 -7 Mq yr _1 as before. No restriction on 7?h 
is necessary, because by the time the gas reaches the inner bound- 
ary, the accretion rate is by definition less than or equal to the Ed- 
dington rate. The results are shown in Figure^ Again, the outburst 
starts in the same way as before, but here, when the large spike in 
surface density and local accretion rate develops at R = 0.15a, a 
large amount of gas is blown away. Unlike in the previous simula- 
tions, this gas never reaches the inner boundary, and although vari- 
ability in the central accretion rate is produced, there is less mass 
involved and the appearance of the outburst profile is simplified. 
This is reflected in the depleted appearance of the disc, especially 
near the inner boundary, at times later than t = 10 orbits (Fig- 
ure |8j. The flares appear to be more regular than in the previous 
case, although there is some evolution in the variability. This is in 
response to the large mass-loss rate at the beginning of the simula- 
tion. Figure[5]shows an expanded section of the mass accretion rate 
through the inner boundary in three 1000-day sections. The period- 
icity of the variations decreases gradually, but by t — 1600 days 
the flares have settles into a quasi-steady state, recurring on a time- 
scale of ~ 80 — 100 days. It should be noted that the time-scale 
of these flares is inextricably linked to the amount of mass lost in 
the disc, and this depends on the value of Ro and the accretion 
efficiency. We discuss this in more detail in Section 4 below. In or- 
der to minimise noise, the central accretion rate is calculated over 
intervals of 50 time-steps during the simulation, rather than con- 
sidering the highly discretized case of counting an integer number 
of particles accreted during a single finite time-step. 50 time-steps 
correspond to 1.3 days in Figure|5| hence we are confident that even 
the shorter time-scale variations that we find in the simulations are 
not attributable to noise. 

The bottom panel of Figure^shows the rate at which mass is 
removed from the disc before it reaches the black hole. After the 
initial burst, the mass loss rate settles down and oscillates around a 
mean value of -M ~ 2 x lO _lo M yr~ 1 . 



4 DISCUSSION 

Our results show that mass loss, irradiation and tidal interactions all 
have a profound effect on the observable properties of long period 




Figure 8. Surface density map of the accretion disc in simulation 3 at t=16 
orbits. The colour density scale is the same as for Figure|5] The inner part 
of the disc is clearly depleted where the wind loss is greatest. 
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Figure 9. Expanded view of the centre panel of Figure Q showing mass 
accretion rate through the inner boundary in simulation 3, beginning at orbit 
10. The data is presented here with time in days (not orbital periods as in 
previous figures) and has been divided into three 1000-day slices to aid 
comparison with FigurelTI 



X-ray transients. In particular, the interplay of all of these effects 
in the outer regions of the accretion disc is able to drive long time- 
scale (i.e. months to years) variability is these objects, and is a pos- 
sible origin for some of the extreme variability of GRS 1915+105. 
The physical process that drives these variations occurs in 
three main stages. In the first stage, after the onset of an outburst, 
the mass accretion rate rises and the disc becomes irradiated out to a 
certain radius Rh . If the Eddington limit holds, this radius is limited 
to a maximum value corresponding to accretion at the Eddington 
rate at the inner boundary of the disc. We have shown that the flow 
through a fixed hot/cold boundary is neither steady nor smooth. 
Gas in the unirradiated part of the disc rapidly crosses the critical 
threshold for the disc instability and the remainder of the accretion 
disc now makes the transition to the hot, high viscosity state. In the 
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third and final stage, the interaction of non-axisymmetric structure 
(here, the spiral waves induced by the tidal stresses of the compan- 
ion star on the disc) with the unsteady hot/cold boundary produces 
a rather inhomogeneous flow. The local variations in density and 
mass accretion rate produce the variations in accretion rate seen in 
our models. A further layer of complexity in the variations is added 
when local mass loss from the disc (here assumed to be in the form 
of a wind) is taken into account. 

The variable which has the most impact on the variability that 
is predicted by this model is the magnitude of Rh.max, the radius 
to which the disc can be irradiated for accretion at the Eddington 
limit. This is sensitive to the geometry of the disc and the illuminat- 
ing X-ray source and to the accretion efficiency. It is probable that 
self-shadowing and radiative warping of the disc will play some 
part in the behaviour, which we have not considered here. If the 
accretion efficiency is particularly low, Rh.max is correspondingly 
closer to the black hole, perhaps at some radius smaller than we 
can resolve for the inner boundary of our disc in these simulations. 
While this will impact upon the time-scale of the variations caused 
by local mass loss - the rate of mass loss depends on the local Ed- 
dington rate, which increases with radius - any variations caused 
by mass loss at smaller radii will occur on shorter time-scales and 
not affect the gross long-term behaviour that we find in simula- 
tion 3. Furthermore, the inhomogeneity that starts to appear near 
the hot/cold boundary rapidly spreads to the whole disc outside 
Bh,max. The triggering factor for the long time-scale variability is 
that the Rh.max boundary remains fixed, regardless of its magni- 
tude. 

A related issue is the duration of the initial high central accre- 
tion rates, which we find at the beginning of simulations 2 and 3. 
Whether such an event is present in the observed X-ray light curve 
in Figure 1 is debatable: the case rests on whether the initial high 
state that decays over the first 100 days or so of the outburst is part 
of a single, coherent event or not. We can, however, comment upon 
the mechanism that produces the initial high state in the simula- 
tions. It is caused by the accretion of material that starts off inside 
Rh.max, and becomes irradiated. It follows that the duration of this 
initial burst will depend on the viscous time-scale at the position 
of Rh.max, which in turn depends on the factors discussed above. 
Here, our choice of accretion efficiency, r\ and MEdd(-Ro) produce 
a value of -Rh.max that is well outside the inner boundary of the 
simulated disc. If nature conspires to make -Rh.max smaller than 
this estimate, the duration of the initial burst would be shorter. 

The postulate that the local wind-loss mechanism must occur 
at small radii and produce short time-scale variations is indeed sup- 
ported by the results of simulation 3. The mean wind loss rate in the 
simulation from radii greater than the radius of the inner boundary 
(Rwind > 3.75 x 10 n cm) is ~ 2 x lQ-^M^yr" 1 . T his is far less 
than th e estimate of 10" 6 to 10" 7 MQyr -1 made by Kot anTet alJ 
(2000) to explain the hydrogen column density inferred from ab- 
sorption line spectra. It is more likely that this kind of mass loss 
rate is generated at much smaller radii in the disc, where MBdd is 
much lower. 

A fully realistic model for the viscous processes at work in 
these discs requires non-ideal magnet ohydrodynamics (MHD), as 
the equivalent value of a is likely to vary with disc radius, vertical 
displacement from the mid-plane and with time. While for numer- 
ical reasons the values of a which we have used in this work are a 
little higher than are usually assumed for these discs, the discrep- 
ancy is not marked. The entire disc spends almost all of its time in 
the high viscosity state during an outburst. MHD calculations have 
shown that the magneto-rotational instability can sustain a mean 



equivalent viscosity parameter of a ~ 0.4 (Tout & Pringle 1992). 
While we use a — 1 here, the viscosity is dominated by the sound 
speed, scaling as v oc ac^. The slight overestimate in a is easily 
compensated by our rather conservative estimate of sound-speed in 
the hot state (described at the beginning of Section 3). 

We close by commenting on other factors that contribute to 
such incredibly high accretion rates during an outburst. If, during 
a long quiescent period the inner regions of the accretion disc are 
absent through evaporation, the accretion disc may build up a very 
large reservoir of mass. This could lead to extremely high accretion 
rates during a subsequent outburst. In the case of GRS1915+105, 
however, we expect an unusually high mass transfer rate from the 
companion star, and this could be a fundamental factor in the abil- 
ity of such systems to achieve and maintain accretion rates near, 
and perhaps beyond, the Eddington limit. Mass-transfer rates of 
the order ~ 10~ 8 MQyr -1 are only achievable in low-mass X-ray 
binaries that contain significantly evolved companion stars, trans- 
ferring mass on a time-scale governed by nuclear expansion of the 
secondary. It is no surprise, then, that GRS 1915+105, with its long 
orbital period and huge accretion disc is one of the few systems in 
which such behaviour has been observed. 
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